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1.  SUMMARY 


The  research  performed  under  this  contract  was  comprised  of  a  number  of  innovative  research 
and  analysis  programs  to  understand  the  fundamental  solar  and  magnetospheric  processes  and 
mechanisms  that  drive  and  control  the  Earth’s  space  environment  in  order  to  understand,  specify, 
anticipate  and  mitigate  the  effects  of  space  weather  on  Air  Force  systems.  As  part  of  our  efforts, 
we  supported  this  program  with  research  and  instrument  development  to  study  the  origins  and 
interplanetary  (IP)  propagation  of  disturbances  that  affect  the  geosphere  and  to  support  the 
mission  and  analyze  the  data  from  the  AFRL  Solar  Mass  Ejection  Imager  (SMEI)  experiment,  a 
heliospheric  image/detector  that  observes  and  tracks  coronal  mass  ejections  from  the  Sun  to 
Earth.  This  program  includes  support  to  AFRL  to  design  and  develop  the  next  generation 
heliospheric  imager  system.  Boston  College  also  provides  scientific  and  technical  collaboration 
to  AFRL/RVB  to  aid  and  enable  the  development  and  validation  of  the  AE-9  and  AP-9  models, 
the  next-generation  of  trapped  electron  and  proton  models  designed  to  replace  the  current  AE-8 
and  AP-8  models  and  provide  users  with  much-needed  improvements  in  accuracy  and  additional 
capabilities.  This  work  also  addresses  the  development  of  enhanced  versions  of  AFGEOSpace 
and  the  Heliospace  analysis  toolset.  Finally,  special  projects  in  support  of  Earth’s  Space 
environmental  studies  are  addressed  throughout  the  duration  of  the  contract.  These  include  the 
analysis  of  various  space  sensors  and  satellite  missions. 

2,  INTRODUCTION 

This  contract  consisted  of  a  number  of  projects  in  support  of  Space  Weather:  Measurements, 
Models  and  Predictions.  In  the  following,  we  introduce  each  project,  describe  the  analysis  and 
provide  the  main  results  of  our  efforts. 

3.  METHODS,  ASSUMPTIONS,  AND  PROCEDURES 

3,1  Solar  Origin  and  Interplanetary  Propagation  of  Geoeffective  Disturbances 

One  of  the  prime  goals  of  this  contract  has  been  to  understand  the  origins  and  interplanetary  (IP) 
propagation  of  disturbances  that  affect  the  geosphere,  and  to  develop,  support  and  analyze  the 
data  from  the  Solar  Mass  Ejection  Imager  (SMEI)  spacecraft  experiment,  designed  to  provide 
early  warning  and  measurement  of  space  hazards.  The  research  studies  have  involved  collecting 
and  analyzing  satellite  and  ground-based  data  sets  of  pertinent  solar,  IP  and  geomagnetic 
phenomena  to  address  3  general  topics:  (1)  The  solar  eruptive  phenomena  which  lead  to 
sporadic  geomagnetic  storms;  (2)  IP  shocks  and  other  signatures  of  solar  ejecta;  and  (3)  The 
characteristics  of  IP  disturbances  which  produce  geomagnetic  storms.  The  most  important  goal 
of  the  studies  is  to  utilize  data  sets  related  to  solar  ejecta  and  their  IP  manifestations  to  define  the 
origins  and  propagation  characteristics  of  geoeffective  disturbances.  The  work  also  involved  the 
development  and  support  of  the  SMEI  spacecraft  experiment  and  analysis  of  the  data  returned. 
SMEI  was  a  proof-of-principle  experiment  for  providing  an  early  warning  of  1-3  days  of  the 
arrival  at  Earth  of  a  coronal  mass  ejection  (CME),  and  measurements  of  crucial  geoeffective 
parameters  of  the  CME. 
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3.1.1  Coronal  Mass  Ejections  (CMEs) 

Coronal  mass  ejections  consist  of  large  structures  containing  plasma  and  magnetic  fields  that  are 
expelled  from  the  Sun  into  the  heliosphere.  They  are  of  interest  for  both  scientific  and 
technological  reasons;  scientifically  because  they  are  responsible  for  the  removal  of  built-up 
magnetic  energy  and  plasma  from  the  solar  corona,  and  because  technologically  they  are 
responsible  for  the  most  extreme  Space  Weather  effects  at  Earth,  as  well  as  at  other  planets  and 
spacecraft  throughout  the  heliosphere.  The  CME  plasma  is  entrained  on  expanding  magnetic 
fields,  which  commonly  have  the  form  of  helical  field  lines  with  changing  pitch  angles,  i.e.,  a 
flux  rope. 

The  recent  solar  minimum  (2008-2009)  between  solar  cycles  23  and  24  was  unusually 
prolonged,  with  record  numbers  of  sunspot-free  days,  record  low  solar  polar  magnetic  fields,  and 
record  high  levels  of  cosmic  ray  flux.  There  were  broad-ranging  terrestrial  responses  to  this 
inactivity  of  the  Sun.  BC  was  involved  in  the  international  Whole  Heliosphere  Interval  (WHI) 
program  to  study  this  minimum  period  and  compare  it  with  previous  minima  [1,  2,  3,  4].  Most  of 
the  ~60  coronal  mass  ejections  (CMEs)  during  WHI  occurred  during  the  first,  active  half  of  the 
interval.  The  origins  and  early  development  of  the  CMEs  and  their  associated  active  regions 
were  analyzed  [5].  The  overall  evolution  of  the  active  regions  and  CMEs  were  analyzed  in  the 
context  of  the  three  rotations  centered  on  CR  2068. 

The  STEREO  Heliospheric  Imagers  (His)  are  part  of  the  SECCHI  suite  of  imaging  telescopes  on 
each  spacecraft  and  view  the  inner  heliosphere  starting  at  an  elongation  of  4°  from  the  Sun.  The 
HI  telescopes  observe  up  to  a  90°  range  in  position  angle,  centered  on  the  ecliptic  and  viewing 
either  east  (HI- A)  or  west  (HI-B)  of  the  Sun-Earth  line.  Combined  with  the  SECCHI 
coronagraphs,  the  His  provide  a  continuous  view  from  the  Sun  to  elongations  of  '-90°  from  the 
Sun-STEREO  line.  During  the  overlapping  period  of  their  missions,  2007  through  2011,  SMEI 
provided  global  context  observations  for  the  SECCHI  His  as  well  as  a  “third-eye”  view  from 
Earth  orbit.  The  combination  of  the  limited-field  HI  and  all-sky  SMEI  observations  form  a 
unique  data  set  that  have  already  provided  fundamentally  new  insights  into  CMEs.  The  BC 
group  and  collaborators  were  involved  in  many  new  studies  using  both  SMEI  and  SECCHI  data 
that  are  helping  develop  a  deeper  understanding  of  the  structure  of  CMEs  and  how  they 
propagate  in  the  heliosphere  [6].  Results  using  3D  global  density  reconstructions  of  the 
propagation  of  CMEs  with  SMEI  and/or  SECCHI  HI  data  were  published  by  [7,  8,  9,  10,  11,  12, 
13,  14,  15]. 

Although  we  know  the  importance  of  CMEs  in  the  evolution  and  propagation  of  geoeffective 
disturbances,  their  solar  origins  remain  obscure.  CMEs  often  exhibit  a  “three-part”  structure 
consisting  of  a  bright  leading  arc  followed  by  a  darker,  low-density  cavity  and  bright  core  of 
denser  material.  With  the  3-D  views  from  SOHO  and  STEREO  it  is  now  clear  that  most  CMEs 
contain  flux  ropes  and  are  likely  associated  with  “magnetic  clouds”  detected  in-situ  by 
spacecraft.  The  frequency  of  occurrence  of  CMEs  tends  to  track  the  solar  sunspot  cycle  in  both 
phase  and  amplitude,  but  during  the  recent  cycle  24  this  relationship  is  diverging  and  is  an  area 
we  are  studying.  Solar  surface  activity  associated  with  CMEs  includes  the  formation  of  dimming 
regions,  long-lived  loop  arcades,  flaring  active  regions,  large-scale  coronal  waves,  and  fdament 
eruptions  and  in  the  corona  with  large-scale,  closed  structures  such  as  coronal  streamers.  During 
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this  period  we  have  studied  all  these  structures,  and  currently  are  investigating  flares,  dimming 
regions  and  arcades  using  SDO  EVE  spectral  and  AIA  imaging  observations  [16].  These  data 
have  great  potential  for  space  weather  forecasting. 

The  trailing  structures  of  CMEs  provide  important  information  on  magnetic  field  line 
reconnection  and  kinematics  of  CMEs.  The  SOHO  EASCO  coronagraphs  have  revealed 
concave-outward,  flux  rope-like  structures  in  nearly  half  of  all  CMEs,  and  trailing  ray-like 
features  suggestive  of  current  sheets  created  during  reconnection.  BC  personnel  have  been  at  the 
forefront  of  eruptive  current  sheet  studies,  were  part  of  the  ISSI  Workshop  on  Eruptive  Current 
Sheets  [17],  and  have  compiled  a  catalog  of  EASCO  post-CME  current  sheets/rays  and  are 
analyzing  them  [18,  19]. 

3.1.2  Solar  Flares 

A  study  was  performed  involving  the  thermal  evolution  of  solar  flares  using  data  from  the 
Extreme  ultraviolet  Variability  Experiment  (EVE)  on  the  recently  launched  NASA  Solar 
Dynamics  Observatory  (SDO).  EVE  has  the  highest  spectral  and  temporal  coverage  of  any  EUV 
experiment  designed  to  study  space  weather  effects.  About  20  flares  observed  over  a  4-month 
period  are  being  studied.  Their  hot  coronal  emissions  ('-'3  MK)  often  have  large  secondary  peaks 
that  appear  many  minutes  to  hours  after  an  eruptive  compact  flare  event.  In  contrast,  the  cool 
coronal  emissions  (-0.5  -  1.5  MK)  usually  dim  immediately  after  the  flare  onset  and  do  not 
recover  until  after  the  delayed  second  peak.  We  refer  to  this  period  of  1-5  hours  after  the  main 
flare  phase  as  the  late  phase.  Our  analysis  of  these  flare  events  suggests  that  the  late  phase 
involves  hot  coronal  loops  near  the  flaring  region,  but  not  directly  related  to  the  original  flaring 
loop  system  [16]. 

3.1.3  The  Solar  Energetic  Particle  Model  (SEPMOD) 

We  initiated  a  study  on  Porecasting  Arrival  Times  of  Shock-Driven  Solar  Energetic  Particles. 

The  objective  of  this  study  was  to  test  and  validate  the  numerical  SEPMOD  model,  which 
estimates  the  solar  energetic  particle  (SEP)  intensity  profiles  for  10  MeV  to  100  MeV  energies  at 
Earth.  SEPMOD  can  be  used  with  a  heliospheric  solar  wind  model,  such  as  Wang-Sheeley-Arge 
(WSA)-Enlil,  together  with  the  geometric  “cone”  description  of  a  CME.  This  modeling  objective 
supports  the  APRE/RVBXS  effort  of  forecasting  the  arrival  times  of  CME  shock-driven  SEPs  in 
order  to  mitigate  their  potentially  harmful  effects  on  Air  Porce  satellites. 

Dr.  Janet  Euhmann  was  consulted  in-depth  with  about  her  SEPMOD  code.  The  current  version 
of  this  model  uses  output  files  generated  by  the  Wang-Sheeley-Arge  (WSA)-Enlil-cone 
modeling  system.  Dr.  Euhmann  gave  training  sessions  from  how  to  use  her  SEPMOD  model  and 
how  to  process  the  output  files.  At  the  conclusion  of  this  meeting  she  shared  her  numerical  code 
and  the  supporting  code  documentation.  The  SEPMOD  code  has  been  installed  and  are  tested  the 
using  a  number  of  test  case  files  provided  by  Dr.  Euhmann.  In  order  to  rigorously  test  the  model 
with  newer  retrospective  SEP  events,  the  model  developer’s  numerical  version  of  the  WSA-Enlil 
heliospheric  solar  wind  model  that  produces  output  files  required  by  SEPMOD  was  used.  Using 
the  current  state-of-knowledge  of  CME-shock  driven  solar  energetic  particles,  a  compilation  was 
created  of  recent  SEP  events  observed  at  earth. 
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3.1,4  Solar  Mass  Ejection  Imager  (SMEI) 

CMEs  have  been  observed  in  white  light  for  decades  with  coronagraphs  viewing  near  the  Sun 
that  observe  Thomson-scattered  brightness  of  sunlight  scattered  from  heliospheric  electrons. 

Over  the  last  decade  heliospheric  imagers,  also  viewing  in  white  light,  were  developed  and 
launched.  The  Solar  Mass  Ejection  Imager  (SMEI),  initially  funded  by  NASA  and  AERL,  was 
launched  on  the  Coriolis  STP  mission  in  January  2003  and  consisted  of  three  baffled  cameras 
whose  3°  X  60°  fields  of  view  were  aligned  in  their  long  dimension  to  achieve  a  combined 
'-'160°-wide  field  of  view  that  scanned  most  of  the  sky  every  102-minute  orbit.  The  majority  of 
SMEI  data  have  an  angular  resolution  of  -0.2°.  Its  global  sky  coverage  permitted  continuous 
imaging  sequences  and,  thus,  a  view  of  heliospheric  structures  as  they  moved  across  the  sky  over 
a  wide  range  of  solar  elongations.  This  type  of  analysis  of  SMEI  data  has  provided  a  unique  new 
perspective  on  heliospheric  observations  of  solar  disturbances  (e.g.,  [20]). 

Modeling  the  propagation  of  CMEs  necessarily  entails  validation  with  heliospheric  imaging. 
However,  the  Thompson  scattered  light  from  CMEs  is  several  orders  of  magnitude  fainter  than 
the  backgrounds  they  traverse  through  the  heliosphere,  so  careful  processing  is  needed  of  the 
images  to  remove  the  backgrounds  but  leave  the  CME  emission.  BC  has  a  legacy  of  working 
with  heliospheric  imaging  as  well  as  CME  tracking  and  modeling.  AERL  and  BC  were 
responsible  for  the  data  pipeline  and  image  construction  of  the  SMEI  Project,  and  producing 
“quick-look”  maps  for  near-real  time  (NRT)  space  weather  forecasting.  The  University  of 
California  at  San  Diego  (UCSD)  SMEI  pipeline  was  designed  to  maintain  as  much  as  possible  of 
the  original  SMEI  data  frame  resolution  and  photometric  precision.  BC  has  worked  closely  with 
the  UCSD  SMEI  group  in  processing  and  analysis  of  SMEI  data  on  CMEs  [10,  11,  12,  7]. 

The  AERE/BC  group  produced  Elexible  Image  Transport  (EITS)  images  in  Aitoff  and  fisheye 
projections  for  NRT  forecasting  with  1°  spatial  resolution,  and  higher  resolution  maps  for 
archiving  purposes.  The  UCSD  database  now  includes  over  4  terabytes  (TB)  of  individual  data 
frame  images,  and  over  6  TB  of  data  products  available  at:  http://smei.ucsd.edu/.  The  BC  team 
has  used  a  set  of  archived  difference  images  to  identify  and  track  CMEs  in  the  SMEI  data,  and 
these  CMEs  have  been  cataloged  using  an  AERL/BC  processing  technique.  Over  the  8.5  year 
lifetime  of  the  mission,  we  have  identified  >400  events  using  this  technique. 

At  BC  many  analyses  of  SMEI  CME  data  was  preformed,  especially  for  individual  events  and  on 
the  solar-cycle  variation  of  CME  rates.  In  addition,  the  backgrounds  in  SMEI  images  provide  a 
wealth  of  useful  astronomical  data  which  was  exploited  by  the  BC  group.  Such  studies  have 
included  the  detection  of  high  altitude  aurorae  [21],  the  disruption  of  comet  tails  [22],  analysis 
and  modeling  of  the  zodiacal  dust,  and  the  long-term  monitoring  of  bright  stars  [23].  The 
zodiacal  dust  appears  as  scattered  sunlight  but  it  has  mostly  been  observed  in  the  infrared  with 
visible-light  models  lacking.  The  SMEI  data  provide  unique  data  to  study  and  model  the 
zodiacal  light  since  SMEI  continually  observed  it  for  nearly  814  years  and  a  realistic  3-D  model 
of  dust  cloud  is  possible.  The  SMEI  database  also  provides  long-term  monitoring  of  bright  stars. 
Many  of  these  stars  are  used  as  calibrators  for  astronomical  surveys,  yet  their  long  term  stability 
as  such  has  rarely  been  investigated.  The  AERL  network  of  calibration  stars  uses  many  of  these 
bright  stars  as  references  in  the  infrared.  This  is  now  being  extended  into  the  visible  and  the 
SMEI  database  is  being  used  to  determine  any  variability  in  the  visible.  Einally,  SMEI  made 
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many  serendipitous  observations  of  satellites  and  space  debris  that  crossed  the  field  of  view  of 
the  cameras.  Some  studies  of  the  nature  and  origin  of  this  debris  has  been  made  by  the  BC 
group,  most  recently  of  the  origin  of  a  debris  swarm  that  was  seen  near  the  cameras  [24], 

As  solar  activity  increased,  SMEI  began  to  observe  more  intense  CMEs/storms  as  they  traversed 
the  heliosphere.  These  observations  were  measured  and  identified  and  movies  of  SMEI  images 
with  the  IDs  were  generated  in  each  case.  Also,  with  more  accurate  timing  and  positions, 
archived  SMEI  data  were  examined  for  launches  that  were  not  observed  in  the  SMEI  final 
processed  images.  In  these  instances,  the  individual  frame  data  were  examined  in  the  eventuality 
that  the  processing  pipeline  for  the  SMEI  imaging  removed  evidence  of  the  launches  as  an 
artifact.  Eor  the  first  few  cases,  one  additional  launch  was  observed.  A  catalog  of  such  “missed” 
launches  has  been  compiled. 

A  SMEI  measuring  tool  was  created  to  provide  measurements  of  CMEs/storms  as  they  traversed 
the  heliosphere.  Minor  updates  were  incorporated  into  the  software  over  the  lifetime  of  the 
mission.  Also  the  software  was  modified  to  provide  a  more  user-friendly  GUI  for  operators  at  the 
Space  Weather  Eorecast  Eaboratory  (SWEE).  It  provided  the  user  with  an  estimate  of  the  launch 
time  for  the  measured  solar  disturbance.  Though  these  were  calculated  in  the  original  program, 
the  results  are  provided  to  the  user  as  a  summary.  Re-evaluation  of  the  state  of  the  health  and 
responsivity  of  the  cameras  over  the  mission  lifetime  was  made  continually  during  the  mission. 
Software  was  written  to  extract  the  photometry  of  individual  stars  on  each  frame  of  data  by 
fitting  each  star  with  an  interpolated  point  response  function  (PRE).  A  catalog  of  the  5800 
brightest  stars  (down  to  6th  magnitude)  was  used  to  locate  the  stars  on  individual  frames.  All  the 
current  pipeline  corrections  used  for  imaging  were  applied  to  each  frame  (dark  correction,  flat 
field,  spike  removal,  optical  path  correction,  etc). 

During  this  period  we  continued  to  maintain  and  update  the  SMEI  catalog  of  CMEs  and  the 
SMEI  publications  list.  A  summary  of  the  mission  goals  and  lessons  learned  was  also  contributed 
to  by  BC  researchers.  After  the  termination  of  the  SMEI  mission  in  201 1,  a  comprehensive 
paper  describing  its  scientific  discoveries  and  legacy  was  published  [25]. 

3.1.5  Space  Weather  Forecast  Lab  (SWFL) 

BC  team  members  were  also  an  integral  part  of  AFRE’s  Space  Weather  Forecast  Eaboratory 
(SWEE).  SWFL  was  created  in  2007  to  integrate  research  advances  in  several  disciplines  to 
provide  new  or  improved  end-user  capability  to  specify  and  forecast: 

1)  satellite  anomaly  hazards 

2)  impacts  of  the  ionosphere’s  state  and  solar  emissions  on  radio  wave  propagation, 
interference  and  scintillation 

3)  upper  atmospheric  changes  that  affect  satellite  drag. 

BC  team  members  attended  routine  meetings  of  SWFL,  led  work  groups  and  contributed  to 
research  and  discussion  in  support  of  SWFL. 
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One  was  using  only  White  Light  (Heliospherie  Imager)  data  for  foreeasting.  The  goal  was  to 
demonstrate  and  evaluate  the  use  of  SMEI  data  to  direetly  estimate  CME  propagation  speed  and 
to  foreeast  shook  and  CME  arrival  times  at  Earth.  At  the  time  of  the  AFRL-East  move  validation 
was  in  progress  oomparing  predioted  arrival  times  with  El /ACE  indioators  and  with  heliospherie 
imagery  to  assess  aoouracy  of  arrival  time  estimates.  Part  of  this  effort  was  to  validate  the 
Tapp  in-Howard  propagation  model  developed  with  AFRL  resouroes.  The  Tapp  in-Howard  model 
assimilates  heliospherie  image  data,  suoh  as  from  SMEI  or  the  STEREO  HI  imagers,  to  track  and 
forecast  CME  arrivals.  The  model  was  used  in  several  studies  for  validation  for  space  weather 
forecasting  [26,  27].  We  also  developed  a  tracking  tool  for  testing  for  operational  use  at  AFWA. 
Both  the  tracking  tool  and  the  Tappen-Howard  (T-H)  model  were  evaluated  from  2010-2011 
using  SMEI  data  to  predict  the  speed  and  time  of  arrival  (ToA)  of  CMEs  at  1  AU  vs  the  time  that 
the  forecast  was  actually  made  [28].  We  were  also  involved  in  improving  predictions  of 
CME/shock  arrival  times  at  Earth  using  the  HAF-v2  solar  wind  model  and  resulting  storm 
strength.  Another  SWFL  task  was  to  evaluate  the  performance  of  several  models,  such  as  WSA- 
Enlil  and  HAF  that  forecast  solar  wind  and  IMF  characteristics  for  a  representative  period  of  a 
solar  cycle. 

BC  personnel  were  also  involved  in  intercomparisons  of  forecasts  of  space  weather  events,  using 
LASCO,  STEREO  and  other  datasets.  During  2010-2011,  STEREO  Group  members  made  near- 
real  time  event  predictions  for  three  events:  8-12  April  2010,  the  30  July-4  August  2010  series 
of  multiple  events,  and  the  15-18  February  2011  X-class  flare-CME  [29,  30].  The  Group 
predictions  for  these  three  sets  of  events  were  based  on  several  models  and  techniques,  including 
geometric  localization  and  polarization,  HI  “J-plots”,  a  NOAA  empirical  model,  the  HAFv2 
model,  the  WSA+Enlil  model,  the  Shock  Time  of  Arrival  (STOA)  model,  the  Empirical  Shock 
Arrival  model,  and  the  Tappin-Howard  model.  For  all  three  event  periods,  the  prediction 
accuracy  averaged  8.0  hr.  with  a  large  range,  an  improvement  over  the  typical  time-of-arrival 
uncertainties  in  the  past  literature  of  1 1-12  hours.  The  lead  or  warning  time  averaged  ~1  day. 
These  results  showed  that  use  of  heliospheric  imaging  data  improves  forecast  performance,  and 
that  a  mission  with  imagers  dedicated  to  space  weather  will  provide  the  best  performance. 

3,1.6,  Hybrid  Space  Weather  Forecasting  and  Modeling 

BC  collaborated  to  support  AFRE’s  efforts  in  developing  a  Next  Generation  Heliospheric 
Imager,  assisting  Janet  Johnston  of  AFRL  and  world-wide  experts  in  two  workshops  held  at 
Sacramento  Peak  Observatory,  NM  on  developing  missions  such  as  to  the  Sun-Earth  E5 
Eagrangian  point.  Such  missions  would  involve  both  heliophysics  and  space  weather  research 
and  forecasting,  with  an  array  of  instruments  including  a  Next  Generation  Heliospheric  Imager 
[31,32]. 

The  development  of  hybrid  space  weather  forecasting  methods  was  performed  under  this  effort. 
Several  categories  were  developed  based  on  precursor  and  advance  warning  time  frame  (e.g. 
forecast  to  nowcast  to  post-cast).  These  time  frames  can  range  from  days  to  hours.  Various 
satellites/observations  were  tagged  as  being  the  most  relevant  for  these  forecast  regimes.  The 
types  of  observations/  precursors  (.e.g.  coronal  holes,  fdament  eruptions,  sigmoids)  and  the 
probabilities  of  an  event  with  each  of  these  categories  were  refined. 
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As  part  of  the  hybrid  forecasting  method,  a  CME  prediction  algorithm  was  developed  based  on 
climatology  data  from  the  previous  solar  cycle.  This  is  based  both  on  sunspot  number  and  CME 
counts  from  Cycle  23.  The  two  data  sets  are  adjusted  and  scaled  for  the  current  cycle,  since  solar 
activity  is  at  a  reduced  rate  compared  to  the  previous  solar  cycle.  Probabilities  are  calculated  for 
the  occurrence  of  CMEs  on  any  given  day  in  the  solar  cycle  as  well  as  the  trends  of  CME 
occurrence  within  the  cycle.  Recurrence  and  persistence  probabilities  of  these  events  which  will 
give  forecasters  near-term  probabilities  as  well  as  historical  ones  were  investigated.  The  output 
of  one  predictor  is  the  monthly  averaged  CME  rate  in  CMEs/day.  Eurther  adjustments  are 
required  since  the  last  solar  cycle,  minimum  to  minimum,  indicated  a  12.5  year  cycle,  rather  than 
the  canonical  1 1  years. 

Since  BC  contributed  the  current  and  persistence  CME  rate  for  2010  to  the  weekly  SWEL 
forecast  meeting,  we  folded  these  statistics  into  the  above  prediction  algorithm  for  a  more 
accurate  shift/scaling  of  the  previous  cycle’s  numbers.  There  is  a  known  bias  in  the  previous 
cycle’s  count  onward  from  about  2005.  Our  goal  was  to  correct  this  bias  in  our  prediction 
algorithm 

Also,  as  part  of  the  hybrid  method,  we  looked  at  the  persistence  CME  rates  from  east-to-west 
limb  events.  That  is,  determining  the  number  of  events  on  any  given  date  and  determining  the 
number  of  events  on  the  opposite  limb  one-half  a  solar  rotation  (+/-  one  day)  later.  This  will 
give  an  indication  of  persistent  solar  activity.  Currently  the  CD  AW  catalog  data  does  not  give 
the  surface  source  regions  for  these  events,  so  a  proxy  of  position  angle  in  the  catalog  is  being 
used  to  indicate  limb  location  in  lieu  of  the  actual  solar  surface  coordinates.  Eor  this  proxy,  it 
was  determined  that  events  with  position  angles  within  20  degrees  of  solar  equator  during  a 
period  of  solar  minimum  activity  reduce  confusion.  Early  results  showed  a  20%  persistence  rate 
from  limb-to-limb. 

3,2  Space  Particle  Modeling:  The  Next  Generation  Trapped  Electron  and  Proton 
Models 

The  radiation  belts  and  plasma  in  the  Earth’s  magnetosphere  pose  hazards  to  satellite  systems 
that  restrict  design  and  orbit  options  with  a  resultant  impact  on  mission  performance  and  cost. 

Eor  decades,  the  standard  space  environment  specification  used  by  the  engineering  community 
was  provided  by  the  NASA  AE-8  and  AP-8  trapped  radiation  belt  models.  With  well  known 
limitations  on  their  validity,  AERL/RVBX  initiated  collaboration  with  the  National 
Reconnaissance  Office  (NRO),  Aerospace  Corporation  and  Eos  Alamos  National  Eaboratory 
(EANE)  to  develop  the  next  generation  trapped  electron  and  proton  models,  to  be  designated 
AE-9  and  AP-9.  These  models  will  replace  the  AE-8  and  AP-8  models  and  provide  users  with 
much-needed  improvements  in  accuracy  and  additional  capabilities  not  possible  with  the  current 
models. 

During  the  course  of  this  contract,  Boston  College  provided  scientific  and  technical  support  to 
AERE/RVBX  to  enable  the  development  and  validation  of  AE-9  and  AP-9  models.  BC  activities 
in  the  project  involved  processing  data  from  sensors  on  multiple  spacecraft  and  organizing  it  into 
a  model  suitable  for  spacecraft  design.  The  primary  objectives  of  the  effort  were  to  1)  Improve 
the  overall  accuracy  of  the  models;  2)  Provide  indicators  of  the  uncertainty  in  the  model  due  to 
natural  variability  and  instrument  uncertainty;  3)  Cover  a  broad  energy  range  including  hot 
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plasma,  relativistic  electrons  and  highly  energetic  protons;  and  4)  Provide  complete  spatial 
coverage. 

The  following  covers  the  prime  activities  provided  by  BC  team  members  for  this  effort.  BC  staff 
member,  Stuart  Huston,  was  instrumental  in  this  project. 

3.2.1  Proton  Flux  Maps,  Our  support  activities  began  with  the  development  of  procedures  for 
binning  the  proton  fluxes  obtained  from  spectral  inversion  of  detectors  on  several  spacecraft  into 
the  standard  grid  established  for  the  AP-9  model.  This  also  involved  developing  techniques  for 
converting  from  one  coordinate  system  (e.g.,  the  invariant  coordinate  system  used  for  the  model) 
to  another  (e.g.,  the  latitude-radius  system  which  is  a  more  intuitive  mapping  in  physical  space). 

The  binning/mapping  exercise  above  revealed  that  the  coordinate  system  based  on  adiabatic 
invariants  results  in  major  inaccuracies  at  low  altitudes.  The  inaccuracy  is  due  to  the  fact  that  the 
trapped  particle  flux  is  controlled  more  by  the  atmospheric  density  than  by  the  magnetic  field, 
and  the  adiabatic  invariance  is  violated.  Several  new  systems  were  investigated,  based  on  the 
altitude  that  a  particle  reaches  as  it  drifts  around  the  Earth.  Data  from  several  spacecraft  have 
been  organized  in  several  ways  in  order  to  determine  which  systems  work  the  best.  Several 
possible  coordinates  have  shown  great  promise  in  providing  much  improved  model  accuracy  at 
low  altitudes,  as  well  as  providing  accurate  temporal  variations  due  to  the  solar  cycle  and  the 
changing  geomagnetic  field. 

The  results  of  this  work  were  included  in  a  briefing  on  proton  flux  mapping  that  was  presented  at 
the  Critical  Design  Review  for  AE9/AP9  held  at  the  Aerospace  Corporation  in  November  2009. 

3.2.2  Data  Sources,  Processiug  aud  Aualysis,  Key  data  for  the  new  models  include  data 
recorded  by  instruments  on  the  TSX5,  S3-3,  Polar,  Azur  and  TacSatd  spacecraft.  Boston  College 
has  supported  the  acquisition,  processing  and  analysis  of  these  satellite  sensor  data  for 
developing  the  radiation  belt  models,  AE-9  and  AP-9 

One  of  our  data  processing  objectives  was  to  generate  and  maintain  the  Compact  Environmental 
Anomaly  Sensor  Experiment  (CEASE)  science  and  ephemeris  databases.  CEASE  consists  of  a 
small  environmental  hazard  detection  system,  which  is  both  lightweight  and  has  low  power 
requirements.  The  total  package  provides  data  on  radiation  dose,  single  event  upsets,  and 
charging.  A  software  system  has  been  developed  that  will  take  the  prime  data  set,  the  CEASE 
Science  fide,  and  parse  the  data  into  individual  time  tagged  data  base  files.  CEASE  data  provided 
throughout  the  duration  of  this  project  was  processed  and  provided  to  AERL. 

Initial  tasks  also  involved  preparing  data  for  our  use.  This  included  cleaning  the  data  to  reveal 
and  remove  potential  contamination  (e.g.,  counts  from  electrons  in  a  channel  that  is  supposed  to 
be  measuring  protons),  dead-time  issues,  and  other  problems  that  might  make  the  data  less 
useful.  Data  cleaning  was  performed  for  the  TSX5/CEASE  data  set;  this  process  revealed  no 
significant  problems.  Eigure  1  shows  an  example  of  plotting  counts  from  a  proton  channel  (D05 
on  the  y-axis)  against  counts  in  a  channel  that  measures  both  electrons  and  protons  (DOl  on  the 
x-axis).  The  plot  reveals  that  D05  is  responding  to  electrons  in  the  outer  zone,  but  simple  filters 
can  eliminate  data  from  the  region  where  the  contamination  occurs. 
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Figure  1 .  Example  of  cross-plotting  counts  in  two  different  channels  to  reveal  potential 

contamination. 


Processing  of  S3-3  Proton  Data 

The  S3-3  satellite  operated  from  1976  through  1979  and  contained  a  proton  telescope  which 
measured  protons  from  80  to  3000  keV  energy.  This  energy  range  is  very  important  for 
phenomena  such  as  surface  dose  on  spacecraft,  but  is  very  difficult  to  measure.  The  S3-3  data 
set  is  nearly  unique  for  this  energy  range  and  coverage  of  the  inner  zone. 

Processing  included; 

Writing  MATLAB  routines  to  read  the  binary  data  files  and  merge  them  with 
adiabatic  invariant  data  calculated  by  AFRL. 

Performing  “data  cleaning”  to  identify  potential  sources  of  contamination,  remove 
noise  from  the  data,  and  identify  any  other  issues. 

Calculate  spectral  correction  factors.  The  S3-3  proton  telescope  has  broad  energy 
channels,  and  the  effective  geometric  factor  of  each  channel  depends  on  the  shape  of 
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the  flux  spectrum  within  that  channel.  A  method  was  developed  to  determine  a 
correction  factor  by  assuming  a  power-law  spectrum  within  the  channel;  the  power- 
law  exponent  was  estimated  using  measured  fluxes  in  adjacent  channels.  This 
correction  was  applied  to  the  measured  fluxes. 

Mapping  fluxes  in  adiabatic  invariant  coordinates.  Figure  2  contains  an  example  of 
these  maps. 


S3-3  P01  115  keV 


Figure  2.  Fluxes  from  the  S3-3  proton  telescope  mapped  in  adiabatic  invariant  coordinates. 
Processing  S3 -3  Electron  Data 

Electron  data  from  S3 -3  was  processed  in  the  same  way  as  the  proton  data.  One  of  the 
significant  aspects  of  this  processing  was  the  “self-calibration”  of  the  S3-3  Magnetic  Electron 
Spectrometer,  which  was  needed  in  order  to  develop  estimates  of  the  uncertainty  of  the 
measurements.  Eor  many  instruments  used  in  this  project,  similar  measurements  are  available 
for  different  sensors  in  the  same  region  of  space  at  the  same  time;  these  are  referred  to  as 
conjunctions.  When  conjunctions  are  available,  it  is  possible  to  cross-calibrate  the  different 
detectors  and  obtain  error  estimates.  Eor  S3-3,  no  comparable  instruments  were  available  for 
conjunctions,  and  so  a  “self-conjunction”  approach  was  taken.  Successive  passes  through  bins 
with  a  time  difference  of  less  than  one  day  were  treated  as  conjunctions,  and  the  flux  from  one 
pass  was  compared  to  the  flux  during  another  pass.  The  standard  deviation  of  the  difference  in 
log(fiux)  between  successive  passes  was  taken  as  the  error  estimate  for  each  energy  channel  of 
the  detector.  One  problem  arose  because  the  higher  energies  exhibited  a  bi-modal  distribution, 
with  a  high  energy  tail  present  in  some  spectra  and  absent  in  others.  Only  spectra  with  the  high 
energy  tail  were  used  for  the  self-calibration  analysis.  However,  future  analysis  of  the  data  could 
reveal  interesting  physics  about  the  presence  or  absence  of  the  tail.  Eigure  3  illustrates  the  self¬ 
calibration  analysis. 
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Figure  3.  Self-calibration  of  S3-3  electron  channels. 


On  the  left  of  Figure  3  is  an  individual  energy  spectra  within  one  bin  for  several  passes  through 
the  bin.  Note  the  two  different  shapes  above  approximately  800  keV.  On  the  upper  right  of 
Figure  3  is  a  comparison  of  flux  (at  1350  keV)  during  one  pass  versus  flux  during  the  subsequent 
pass.  Note  the  “forks”  at  low  fluxes,  which  are  due  to  the  high  energy  tails  in  the  energy  spectra. 
A  histogram  of  the  difference  in  log  flux  between  successive  passes  is  on  the  lower  right  of  the 
figure.  The  standard  deviation  is  indicated  by  the  dashed  red  lines,  and  this  is  taken  as  an 
estimate  of  the  uncertainty  in  the  measurements. 


Due  to  its  orbit  and  instrumentation,  the  S3-3  spacecraft  provided  a  unique  data  set  for  inner- 
zone  electrons.  We  had  hoped  to  use  this  data  in  Version  1.0  of  AE9.  However,  comparisons  of 
S3-3  data  with  AES  and  other  data  sets  in  AE9  showed  that  fluxes  measured  by  S3-3  were 
significantly  higher  than  the  other  data.  We  therefore  investigated  why  S3-3  might  be  so 
different  from  the  other  data. 


Part  of  this  investigation  involved  a  paper  by  [33]  which  contained  data  from  S3-3.  In  particular, 
we  noted  a  figure  which  compared  data  from  S3 -3  with  other  data  sets  and  with  AES  predictions. 
The  data  in  this  figure  were  processed  by  A.  L.  Vampola,  the  Principal  Investigator  (PI)  on  the 
S3-3  particle  detector.  The  data  that  we  processed  for  AE9  came  from  Vampola,  and  we  used 
software  provided  by  him  to  process  the  data.  However,  in  comparing  the  fluxes  from  our 
process  to  those  from  the  1994  paper  [33],  our  fluxes  are  a  factor  of  10  -  20  higher  than  those  in 
the  paper.  Figure  4  illustrates  these  comparisons  where  the  AES  predictions  and  S3-3  data  are 
shown  in  black  [33];  the  S3-3  data  in  red  are  from  our  data  analysis  (using  s33_read)  and  S3-3 
data  in  green  are  s33_read  fluxes  divided  by  16.  A  review  of  the  data  processing  revealed  some 
small  changes  between  1994  and  the  current  software.  These  changes  explain  slight  differences 
in  the  spectral  shape,  but  not  the  factor  of  10  -  20.  Because  we  could  not  explain  this 
discrepancy,  we  decided  to  omit  the  S3-3  data  from  AE9. 
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Figure  IV*1 


Figure  4.  Comparison  of  S3-3  electron  data. 


Processing  Polar  IPS  Proton  Data 

The  Imaging  Proton  Sensor  (IPS)  on  the  Polar  spacecraft  provides  data  on  proton  fluxes  from 
approximately  0.01  MeV  to  2.0  MeV.  The  highly  elliptical  polar  orbit  covers  almost  the  entire 
trapped  radiation  region.  Processing  included; 

Writing  MATLAB  routines  to  read  the  binary  data  files  and  merge  them  with 
adiabatic  invariant  data  calculated  by  AFRL. 

Performing  “data  cleaning”  to  identify  potential  sources  of  contamination,  remove 
noise  from  the  data,  and  identify  any  other  issues. 

Mapping  fluxes  in  adiabatic  invariant  coordinates.  Figure  5  contains  an  example  of 
these  maps. 

A  major  part  of  the  Polar  data  processing  involved  cross-calibration  with  other  detectors. 

The  cross-calibration  was  done  with  the  EPAM  instrument  on  the  ACE  spacecraft  and  the 
MICS  instrument  which  was  also  on  the  Polar  spacecraft.  These  two  instruments  overlapped 
with  the  upper  and  lower  energy  ranges  of  IPS,  respectively.  The  cross-calibration  process 
revealed  some  contamination  issues  where  higher  energy  protons  were  contaminating  the 
lower  energy  proton  measurements.  After  filtering  out  these  contaminated  measurements, 
the  IPS  data  agreed  well  with  both  MICS  and  EPAM  over  the  energy  range  0.05  -  2.0  MeV. 
Eigure  6  shows  a  comparison  between  IPS  and  EPAM.  The  cross-calibration  also  allowed  us 
to  estimate  measurement  errors  in  the  IPS  data. 
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Figure  5.  Cross  plot  of  Polar/IPS  fluxes  vs.  ACE/EPAM  fluxes  at  six  different  energies. 


3,2,3  Spectral  Inversion  Techniques,  Our  previous  speetral  inversions,  necessary  to  derive 
differential  flux  from  detectors  with  wide  energy  response,  used  an  assumed  spectral  function 
with  a  power-law  spectrum  from  10-100  MeV  and  an  exponential  tail  above  100  MeV.  This 
spectral  shape  was  based  on  the  APS  model,  but  it  involved  parameters  which  were  somewhat 
arbitrary.  Another  approach,  based  on  Principal  Component  Analysis  (PCA)  proved  useful  in 
inverting  electron  spectra,  so  a  similar  approach  was  developed  for  protons.  The  approach  was 
based  on  an  analytical  model  by  Selesnick  [34]  and  covers  energies  from  10-1000  MeV. 

Principal  components  were  derived  from  Selesnick’ s  model  and  were  used  to  obtain  spectra  from 
the  TSX5/CEASE  data  set.  Eigure  6  shows  an  example  of  spectra  obtained  from  the  PCA  model 
(blue)  and  the  previous  analytical  model  (red).  Dashed  lines  indicate  confidence  limits. 
Inversions  based  on  the  PCA  approach  give  slightly  better  agreement  with  the  data,  and  this 
approach  was  used  to  re-do  the  spectral  inversions  for  the  TSX5/CEASE  data  set,  as  well  as 
other  data,  for  the  Version  1.0  release  of  AP9. 


Eigure  6.  Comparison  of  spectral  inversions  of  TSX5/CEASE  using  PCA  approach  and 

analytical  approach. 
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Further  improvements  were  made  to  the  PCA  model  to  extend  the  energies  to  as  low  as  5  MeV. 
These  improvements  were  made  to  enable  the  model  to  be  used  for  speetral  inversion  of  data 
from  other  instruments  on  the  HEO  spaeeeraft.  MATLAB  seripts  were  developed  to  use  the 
PCA  model  to  perform  speetral  inversion  of  this  data.  These  seripts  were  transferred  to  AFRL 
personnel  to  perform  the  aetual  inversions. 

Combined  Spectral/ Angular  Inversion. 

The  flux  maps  produeed  for  AE9/AP9  map,  the  flux  of  loeally-mirroring  partieles  at  a  given 
loeation.  Many  deteetors  (ineluding  TSX5/CEASE)  have  a  wide  aceeptanee  angle  and  thus 
measure  a  wide  range  of  partiele  piteh  angles.  Therefore,  a  method  is  required  to  estimate 790 
from  the  wide-angle  response.  Our  previous  inversions  used  an  after-the-fact  approach  in  which 
the  energy  inversion  was  performed  first,  then  a  correction  factor  was  calculated  by  integrating 
an  assumed  particle  angular  distribution  with  the  known  angular  response  of  the  detector.  This 
approach  was  an  approximation,  but  it  had  the  advantage  that  we  were  able  to  determine  where 
large  corrections  were  necessary. 

In  preparing  to  perform  new  inversions  for  Version  1  of  the  model,  we  developed  an  approach 
which  calculated  an  effective  response  function  for  each  detector  channel  calculated  by 
integrating  the  assumed  particle  angular  distribution  with  the  known  angular  response  of  the 
detector.  This  integration  was  done  before  the  spectral  inversion  and  in  principle  should  result 
in  more  accurate  fluxes.  In  reality,  however,  this  process  resulted  in  many  instances  where  no 
inversion  was  obtained  because  the  effective  response  function  significantly  changed  the  energy 
range  of  the  detector.  It  was  therefore  decided  not  to  use  this  approach  for  the  next  round  of 
inversions. 

We  re-visited  the  approach  after  the  release  of  Version  1.0;  at  that  time  we  were  able  to  develop 
improved  pitch  angle  distribution  models  and  modify  the  inversion  routines  to  make  better  use  of 
the  effective  response  functions.  At  that  time,  progress  was  made  towards  developing  a 
procedure  to  perform  combined  spectral  and  angular  inversions  of  wide-angle  detector  data. 
Proton  pitch  angle  distributions  (PADs)  were  developed  based  on  the  AP-9  model  and  MATLAB 
functions  were  developed  to  integrate  the  ambient  PAD  with  the  angular  response  of  the  detector 
to  obtain  an  effective  response  function  which  can  be  used  in  the  inversion  procedure.  Several 
integration  procedures  were  investigated;  eventually  a  very  simple  procedure  using  the 
rectangular  rule  was  selected,  primarily  because  speed  is  an  essential  feature.  The  integration 
functions  still  need  to  be  tested  and  integrated  with  the  rest  of  the  inversion  procedures. 

3.2,4  Development  of  Flux  Templates,  In  populating  the  run-time  data  tables  for  the  AE- 
9/ AP-9  models,  individual  data  sets  are  combined  into  one  data  base.  None  of  these  individual 
data  sets  covers  the  full  spectral  or  spatial  extent  of  the  final  model.  In  order  to  perform  better 
interpolation  and  extrapolation  of  the  data  sets,  we  are  developing  a  series  of  flux  templates 
which  specify  how  we  expect  the  flux  to  vary  as  a  function  of  energy  and  space. 

Elux  templates  are  being  developed  based  on  a  combination  of  spacecraft.  One  of  the  challenges 
is  extrapolating  to  regions  where  data  do  not  exist.  This  has  been  accomplished  by  fitting  a  sum 
of  exponential  functions  to  the  energy  spectrum  in  each  bin.  The  exponentials  allow  reasonable 
extrapolation  in  the  energy  domain.  The  fitting  parameters  are  then  smoothed  in  the  spatial 


14 

Approved  for  public  release;  distribution  is  unlimited. 


domain,  allowing  a  more  reasonable  spatial  extrapolation.  The  smoothed  parameters  are  then 
used  to  generate  flux  maps  at  the  AP-9  energies;  these  constitute  the  templates.  Figure  7  shows 
an  example  of  a  template.  Dots  indicate  where  data  were  available. 


Figure  7.  Template  mapping  proton  flux  at  1  MeV  based  on  a  composite  of  3  spacecraft. 

3,2,5  TacSat4,  TacSatd  (Tactical  Satellite  4)  is  a  small  spacecraft  launched  by  NRL  in 
September  201 1;  it  carries  a  CEASE  instrument  similar  to  that  on  TSX5.  The  data  from  this 
instrument  will  be  very  useful  in  updating  the  flux  maps  in  AE9/AP9.  We  are  working  with  Dr. 
Chad  Lindstrom  at  AFRE  in  performing  spectral  inversion  on  the  CEASE  data.  One  issue  is  that 
in  order  to  extend  the  inversions  to  lower  energies  (down  to  about  1  MeV),  we  are  including 
several  channels  with  complicated  response  functions.  Eigures  8  and  9  show  the  response  of  one 
of  these  channels  to  both  ‘hard’  and  ‘soft’  energy  spectra.  Red  curve  (left  axis)  shows  the  input 
proton  spectrum.  Blue  curve  (right  axis)  shows  the  channel  response  function  for  isotropic 
fluxes.  Green  bar  indicates  the  energy  range  representing  80%  of  the  total  counts;  vertical  tick 
shows  the  effective  midpoint  energy.  It  can  be  seen  that  the  main  response  of  these  channels  is 
highly  dependent  on  the  spectral  shape. 


Eigure  8.  Response  of  one  TacSatd  CEASE  channel  to  a  ‘hard’  proton  spectrum  (spectral  index 

-1). 
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Figure  9.  Response  of  one  TaeSat4  CEASE  ehannel  to  a  ‘soft’  proton  speetrum  (speetral  index 

-4). 


Additional  work  continued  using  TacSat-4  data.  The  following  summarizes  these  efforts: 

•  One  issue  is  that  in  order  to  extend  the  inversions  to  lower  energies  (down  to  about  1  MeV), 
we  are  including  several  channels  with  complicated  response  functions.  These  channels 
respond  to  both  low-energy  protons  and  high-energy  protons.  We  developed  a  criterion 
based  on  count  rates  in  two  different  channels  to  tell  us  when  the  detector  was  indeed 
responding  to  low-energy  (<  10  MeV)  protons.  Eigure  10  illustrates  this  criterion;  in  general, 
CEASE  responds  to  <  10  MeV  protons  beyond  ~  1.8,  i.e.,  in  the  “slot”  region.  This 
criterion  will  greatly  improve  our  ability  to  interpret  the  CEASE  data. 

•  TacSat-4  Flux  Comparisons.  In  January  and  March  2012,  a  pair  of  moderate  solar  particle 
events  (SPEs)  occurred.  These  events  are  useful  for  the  calibration  of  proton  detectors 
because  they  bathe  the  outer  magnetosphere  (beyond  Lm  ~  6.6)  in  a  relatively  uniform  flux  of 
medium-  to  high-energy  protons.  The  GOES  spacecraft  in  geosynchronous  orbit  are  an 
accepted  standard  for  measuring  these  SPEs;  detectors  on  TSX5,  HEO,  and  other  non- 
geosynchronous  spacecraft  have  been  cross-calibrated  against  GOES  in  the  past  by  studying 
“conjunctions”  when  the  individual  spacecraft  are  at  high  altitude  and/or  high  latitude.  We 
performed  a  similar  analysis  for  TacSat-4/CEASE.  Figure  1 1  shows  energy  spectra  from 
GOES  and  CEASE;  the  agreement  is  quite  good.  This  result  gives  us  confidence  in  the 
absolute  magnitudes  of  the  proton  fluxes  derived  from  CEASE  measurements. 

•  TacSat-4  Attitude  Determination.  The  proton  fluxes  in  Earth  orbit  are  typically  very 
anisotropic;  for  most  detectors  it  is  therefore  critical  to  know  the  look  angle  of  the  sensor 
relative  to  the  magnetic  field.  For  TacSat-4,  the  CEASE  detector  looks  in  the  -X  direction  in 
spacecraft  body  coordinates.  In  order  to  determine  the  look  direction  in  magnetic 
coordinates,  it  is  necessary  to  use  quaternions  supplied  in  the  spacecraft  attitude  control 
system  (ACS)  data  files;  these  calculations  provide  the  look  direction  in  inertial  coordinates. 
This  vector  can  then  be  combined  with  a  magnetic  field  vector  determined  by  a  magnetic 
field  model  to  determine  the  pitch  angle  of  the  detector  Foresight.  We  developed  MATLAB 
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routines  to  read  the  ACS  data  files,  perform  the  neeessary  calculations,  and  merge  the  pitch 
angle  data  with  the  CEASE  flux  data. 

•  TacSat-4  Pitch  Angle  Distributions.  TacSat-4’s  orbit  takes  it  through  essentially  the  same 
regions  of  space  several  times  per  day.  Over  the  first  several  months  of  data  the  pitch  angle 
in  a  given  spatial  “bin”  has  varied  over  a  range  of  about  10  to  90  degrees.  It  has  therefore 
been  possible  to  generate  nearly  complete  pitch  angle  distributions  (PADs).  Figure  12  shows 
a  set  of  PADs  at  several  energies  in  one  bin  near  the  magnetic  equator.  These  PADs  can  then 
be  integrated  to  determine  the  omnidirectional  flux  for  direct  comparison  to  AP9. 


Figure  10.  Plot  of  count  rates  in  CEASE  EB(5,0)  vs.  LB-SUM3. 

When  LB(5,0)  is  more  than  2  x  EB-SUM3  (indicated  by  upper  diagonal  line),  CEASE  channels 
are  responding  to  <  10  MeV  protons.  Color  scale  indicates  Mcllwain  E  value  of  measurements. 


Figure  1 1 .  Comparison  of  CEASE  proton  fluxes  (blue)  with  fluxes  measured  by  GOES- 13 

during  SPE  in  January  2012. 
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Red  squares  indicate  uncorrected  GOES  fluxes;  red  circles/line  indicate  GOES  flux  corrected  for 
spectral  response. 
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Eigure  12.  Pitch  angle  distributions  obtained  from  CEASE  Oct  2011  through  May  2012. 

•  TacSat-4  omnidirectional  daily  fluences.  Data  from  the  CEASE  instrument  on  TacSat-4  has 
been  used  to  model  the  response  of  solar  arrays  to  the  space  environment.  In  order  to  drive 
the  solar  array  response  models,  it  was  necessary  to  determine  daily  omnidirectional  proton 
fluences  from  the  CEASE  data. 

The  determination  of  daily  omnidirectional  fluences  was  complicated  by  several  factors: 

-  The  spectral  inversion  technique  used  on  CEASE  determines  a  unidirectional  flux  which 
reflects  the  average  over  the  CEASE  opening  angle  and  near  the  pitch  angle  of  particles 
entering  along  the  detector  boresight.  Some  way  is  needed  to  estimate  the  full 
omnidirectional  flux  from  these  measurements. 

-  There  are  gaps  in  the  telemetered  data;  these  gaps  may  last  from  a  few  minutes  to  a  large 
fraction  of  a  day.  On  days  where  there  are  significant  gaps,  it  is  impossible  to  determine 
a  true  daily  fluence,  and  a  method  for  filling  the  gaps  is  required. 

Given  these  factors,  it  would  not  be  possible  to  determine  accurate  omni-directional  fluxes  or 
fluences  from  the  CEASE  data  directly.  Eortunately,  there  are  other  factors  working  in  our 
favor: 

-  Over  the  course  of  six  months  or  so,  the  proton  fluxes  in  most  locations  remained 
relatively  stable. 

-  Over  this  time  period,  due  to  orbital  mechanics  and  the  pointing  scheme  for  TacSat-4,  in 
most  bins  the  CEASE  boresight  covered  a  range  of  pitch  angles  sufficient  to  obtain  an 
approximate  pitch  angle  distribution. 

A  procedure  was  developed  which  fills  in  gaps  in  the  data  with  the  average  omni-directional  flux 
in  the  appropriate  bins.  The  primary  assumption  involved  is  that  the  pitch  angle  distribution  in  a 
given  bin  on  a  given  day  matches  the  average  pitch  angle  distribution  (averaged  over  the  entire 
period  of  study)  in  that  bin.  This  procedure  was  used  to  determine  daily  fluences  for  the  first  six 
months  of  the  mission.  Eigure  13  shows  the  average  daily  flux  of  50  MeV  protons  mapped  in  K- 
O  space  for  one  day  of  the  mission.  Superimposed  on  the  flux  map  is  the  spacecraft  ephemeris. 
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The  color  code  indicates  the  daily  average  flux  of  50  MeV  protons.  The  black  lines  show  the 
TacSat-4  ephemeris  for  that  day.  Note  the  high  fluxes  for  logioO  -0.3.  These  are  due  to  an 
ongoing  solar  proton  event. 
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Figure  13.  Daily  average  fluxes  for  24  January 
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•  TacSat-4  Data  Analysis.  In  addition  to  the  work  described  above,  some  more  detailed 
analysis  of  the  TacSat-4/CEASE  data  was  performed.  Eigure  14  compares  energy  spectra 
measured  by  TacSat-4  with  AP9  model  predictions  and  some  of  the  raw  data  which  went  into 
AP9.  It  can  be  seen  that  TacSat-4  measurements  agree  reasonably  well  with  AP9  for 
energies  above  about  20  MeV  and  below  2  MeV.  In  between,  however,  TacSat-4 
measurements  are  considerably  higher  than  AP9  or  the  other  data  sets.  Eor  example,  at  4.3 
MeV,  TacSat-4  is  about  15  times  higher  than  CRRES  and  6  times  higher  than  AP9. 

Figure  15  shows  L-profiles  (at  the  magnetic  equator)  for  CRRES  data,  the  AP9  model,  and 
TacSat-4.  Again,  TacSat-4  is  considerably  higher  than  CRRES  or  AP9. 

Figure  16  shows  time  histories  of  the  count  rate  in  several  CEASE  channels  which  respond  to 
<10  MeV  protons.  Also  shown  is  a  fit  to  the  trend  for  the  LB-B32  channel,  which  responds 
to  2.5  -  4  MeV  protons.  There  is  a  definite  decrease  in  the  count  rate  (and  hence  the  flux) 
over  time.  Based  on  the  inverted  flux  data,  the  flux  Of  4.3  MeV  protons  decreased  by  about 
33%  over  the  course  of  six  months.  Extrapolating  this  rate  of  decrease,  the  TacSat-4  flux 
would  reach  the  levels  measured  by  CRRES  in  about  3.5  years.  Thus,  it  may  be  that  TacSat- 
4  is  measuring  a  long-lived  transient  event.  Continued  monitoring  by  TacSat-4,  combined 
with  measurements  from  additional  spacecraft,  such  as  TWINS  and  RBSP,  would  be  very 
useful  in  explaining  this  phenomenon. 
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Figure  14.  Energy  speetra  from  Polar/IPS,  CRRES/Protel,  TaeSat4,  and  AP9  in  one  bin  near  the 

magnetic  equator. 


Figure  15.  L-profiles  of  10.7  MeV  flux  measured  by  CRRES  and  TacSat4,  compared  with  AP9 

predictions. 
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Figure  16.  Time  history  of  count  rates  in  several  TacSat4/CEASE  low  energy  channels. 
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AP9  Pitch  Angle  Distributions.  For  our  next  release  of  AE9/AP9,  we  are  planning  to  include 
data  from  TacSat-4/CEASE  and  to  re-process  data  from  TSX5/CEASE  using  a  combined 
angular/energy  inversion  technique  which  includes  both  the  angular  response  and  the  energy 
response  of  the  instrument.  In  order  to  do  so,  we  need  a  model  of  the  pitch  angle  distribution  of 
the  ambient  proton  flux.  We  have  used  the  directional  flux  capability  of  AP9  to  derive  pitch 
angle  distributions  as  a  function  of  Mcllwain  L  parameter  and  energy.  Eigure  17  shows  an 
example  of  the  pitch  angle  distributions  at  a  number  of  energies  at  1^=1. 5.  Here,  y=sinaeq, 
where  Oeq  is  the  equatorial  pitch  angle,  and  yLc  is  the  value  of  y  corresponding  to  the  loss  cone. 
These  pitch  angle  distributions  have  been  tabulated  and  will  be  included  in  updated  inversion 
procedures  for  processing  future  data  from  CEASE  and  other  detectors. 


Eigure  17.  Pitch  angle  distributions  derived  from  AP9  for  Em=l-5. 

3,2,6  Solar  Particle  Event  Database,  In  analyzing  data  for  AE9/AP9,  we  have  attempted  to 
remove  data  taken  during  solar  particle  events  (SPEs),  because  the  high  energy  protons  present 
during  these  events  can  contaminate  measurements.  We  have  been  using  two  different  event 
lists  based  on  GOES  data,  one  for  >100  MeV  protons  and  one  for  >  10  MeV  protons.  These 
lists  contained  events  from  1986  through  2006.  Some  of  the  AE9/AP9  data,  however,  were 
taken  before  1986,  and  new  data  (e.g.,  the  TacSat-4  data)  are  being  taken.  In  order  to  develop  a 
relatively  consistent  list  of  events,  a  new  list  has  been  compiled.  This  list  is  based  primarily  on 
the  SEPEM  data  set  (http://dev.sepem.oma.be/)  which  contains  data  from  1973  through  2006. 
The  SEPEM  list  is  slightly  different  in  that  it  uses  the  differential  flux  in  the  7.23  -  10.46  MeV 
channel,  while  the  previous  list  used  integral  fluxes  in  a  >  10  MeV  or  >  100  MeV  channel. 
Comparisons  of  the  two  lists  showed  that  almost  all  the  events  in  one  list  were  contained  in  the 
other;  timings  of  the  events  were  slightly  different.  The  SEPEM  list  was  updated  with  data  from 
GOES  for  2007  -  2012.  A  MATEAB  data  set  was  created  and  a  MATEAB  function  was  written 
to  return  a  flag  indicating  whether  or  not  a  SPE  was  in  progress  at  a  given  time. 
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3.2,7  Relativisitic  Electron  Prediction  Models  Comparisons,  Extended  periods  of  relativistie 
eleetron  intensity  at  geosynehronous  orbit  ean  ereate  severe  deepeharging  hazards  for  satellites. 
Over  the  last  20  years  a  number  of  models  have  been  developed  to  prediet  eleetron  flux  levels 
using  solar  wind  and  geomagnetie  indiees  as  inputs.  We  analyze  the  results  of  several  of  these 
ineluding  the  Relativistie  Eleetron  Eoreeast  Model  based  on  the  linear  predietion  filter  teehnique, 
a  neural  network  algorithm,  and  the  physies  based  diffusion  method.  Analyses  using  the  methods 
of  simple  persistenee  and  recurrenee  (based  on  the  27  day  solar  rotation)  are  also  included  as 
performance  baselines.  Comparisons  are  made  to  the  GOES>  2MeVelectron  flux  to  determine 
which  model  or  method  gives  the  best  +1,  +2,  and  +3  day  forecasts  of  average  daily  flux  during 
the  interval  1996-2008.  Model  inputs  include  combinations  of  SKp,  the  daily  average  solar  wind 
speed,  and  daily  average  >  2  MeV  electron  fluxes  for  one  day  or  multiple  days  prior  to  the 
forecast  days  of  interest.  Prediction  efficiencies  are  calculated  for  6  month  intervals.  After 
evaluating  all  the  models,  there  was  no  clear  winner;  each  model  did  well  at  different  phases  of 
the  solar  cycle.  All  models  perform  their  best  during  the  inclining  phase  of  solar  minimum  but 
not  as  well  during  solar  maximum  and  the  declining  phase  of  solar  minimum.  While  persistence 
is  respectable  for  +1  day  prediction,  models  clearly  give  superior  +2  and  +3  day  predictions  and 
should  be  used  to  obtain  those  forecasts.  The  results  of  this  study  were  published  is  the  Space 
Weather  Journal  [35]. 

3,3  AFGEOSpace  Development 

One  of  our  goals  under  this  contract  was  to  contribute  to  the  maintenance  of  AEGEOSpace2.5 
and  to  the  new  version  AEGEOSpace3.0. 

AE-GEOSpace  is  a  user-friendly,  graphics-intensive  software  program  bringing  together  many  of 
the  space  environment  models,  applications,  and  data  visualization  products  developed  by  the 
Air  Eorce  Research  Eaboratory  and  others  in  the  space  weather  community.  AE-GEOSPACE 
currently  serves  as  a  platform  for  rapid  prototyping  of  operational  products,  scientific  model 
validation,  environment  specification  for  spacecraft  design,  mission  planning,  frequency  and 
antenna  management  for  HE  communications,  and  post-event  anomaly  resolution. 

AE-GEOSpace  provides  common  input  data  sets,  application  modules,  and  1-D,  2-D,  and  3-D 
visualization  tools  to  all  of  its  models.  Available  graphical  tools  include  animation,  annotation, 
axes,  coordinate  probes,  coordinate  slices,  detector  cones  (e.g.,  from  satellites).  Earth,  emitters 
(e.g.,  radar  fans  and  communication  domes),  field-lines  (e.g.,  geomagnetic),  global  inputs  (e.g., 
Kp  and  Dst  indices),  grids,  isosurfaces,  links  (e.g.,  ground-to-satellite),  orbit  probe  (data  along 
orbit  tracks),  orbit  slice  (data  in  orbit  plane),  plane  slice  (data  in  arbitrary  plane),  ray  trace 
(through  ionosphere  electron  densities),  satellites,  stars,  stations,  and  volume  (global  3-D 
rendering  of  data  sets).  The  software's  "dynamic"  mode  enables  the  simultaneous  display  of 
output  from  multiple  environment  models  as  a  function  of  time  over  user-specified  intervals. 

Boston  College  efforts  to  support  AEGEOSpace  2.5  primarily  included  programming,  bug 
fixing,  graphics  development,  resolving  differences  on  program  results  between  operating 
systems,  library  testing  and  responding  to  AERL  requests  for  clarification  and  further 
developments.  One  of  the  problems  with  this  and  earlier  versions  of  AEGEOSpace  were  that 
each  new  component  model  and  its  accompanying  GUI  needed  to  be  compiled  directly  into 
AEGEOSpace.  This  required  the  entire  AEGEOSpace  program  to  be  updated  when  a  single 
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component  model  was  added.  Similarly,  when  models  or  their  interfaee  were  revised,  the  entire 
AFGEOSpaee  program  had  to  be  updated.  During  the  eourse  of  this  eontract,  BC  expended  mueh 
support  to  improve  these  issues  in  the  development  of  AFGEOSpaee  Version  3.0.  These  efforts 
ineluded  eontributing  to  the  development  of  a  modular  approaeh  with  plug-ins  ereated  using 
dynamie  linked  libraries  for  Windows  users  or  shared  libraries  for  Einux  and  OSX  users.  These 
ehanges  allowed  AFGEOSpaee  3.0  to  dynamieally  load  the  module  at  runtime  and  when  a 
eomponent  model  is  ehanged,  AFGEOSpaee  needs  only  to  be  restarted  to  observe  the  changes. 
Sinee  the  models  are  loaded  at  runtime  only  those  models  of  interest  to  the  user  are  required. 
AFGEOSpaee  Version  3  also  ineorporated  a  platform  independent  graphieal  user  interfaee 
(GUI).  That  is  ereated  in  an  external  program  and  then  stored  as  an  XME  file.  This  XML  file  is 
then  loaded  by  AFGEOSpaee  when  the  plug-in  module  is  loaded.  Sinee  the  XML  file  is  a  plain 
text  deseription  of  the  GUI,  it  is  platform  independent  and  only  one  XML  file  is  needed  to  work 
on  Windows,  Linux  or  OSX. 

All  of  our  efforts  for  AFGEOSpaee  2.5  and  3.0  were  coordinated  with  other  the  AFRL  program 
manager  for  AFGEOSpaee  and  with  other  AFRL  eontraetors. 

Some  of  the  work  addressed  by  Boston  College  ineluded: 

•  Development  on  an  interfaee  builder.  This  program  allows  programmers  to  interaetively 
build  GUI's  for  their  plug-ins  for  AFGEOSpaee.  A  "plug-in"  is  a  eode  library  that  is 
dynamieally  loaded  into  AFGEOSpaee  at  runtime. 

•  Implementing  window  ereation  and  ereating  a  menu  bars 

•  Updating  the  GUI  for  the  SNRTACS  model. 

•  GuiBuilder  (  program  that  allows  user  to  build  interfaees  ) 

•  Fixed  issues  when  loading  a  previously  saved  file  objects  would  not  appear  correetly. 

•  Various  debugging  to  the  driver  for  the  SNRTACS  model 

•  SNRTACS  model  diserepaneies  when  run  on  Windows  and  Linux  operating  systems 

•  Fixes  to  the  EPHEMERIS  model  and  eorresponding  graphieal  display 

•  METEOR  IMPACT :  Open  Model  not  working 

•  METEOR  SKY  MAP:  Save/Open  Model  loads  GUI  settings  and  the  eorreet  data  for 
display,  but  any  entries  in  the  "Storms"  box  are  not  retained  -  thus  rerunning  the  loaded 
module  does  not  include  the  missing  storm. 

•  METEOR  IMPACT-APP:  Open  Model  not  working 

•  METEOR  IMPACT-APP:  Freezes  with  valid  Optional  Storm  added  in  fresh  AF¬ 
GEOSpaee  session.  Results  in  "invalid  storm  peak"  message. 

•  Raytraee  App:  Save  Model  message  “AppRayTraeel_Frequeney_time_02  Def  Dim:0 
NC  UNLIMITED  size  already  in  use 

•  Error  Message:  Model  Save  19:  NC-UNLIMITED  size  already  in  use” 

•  Open  Model  message:  Error  Message:  Variable  not  found 

•  Viewport  Menu:  Projeetion:  Attempting  to  use  the  Projection  option  to  change  the 
dimension  of  speeial  1-D  plot  windows  ereated  by  LET-APP,  RAYTRACE-APP,  and 
PPS  causes  erash. 

•  LET-APP:  Use  of  "Trapped  Protons:  CRRESPRO  Quiet"  or  "Active"  results  in  no  output 
ehange  and  error  message. 
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•  Crresele-App:  Exceeded  NETCDF  maximum  number  of  dimensions. 

•  Changes  to  the  management  of  the  dynamie  loading/unloading  of  libraries,  (i.e. 
‘plugins’) 

•  Bug  fixes  to  the  display  of  color  plots 

•  Generating  a  more  unified  build  environment 

•  Heliospaee  changes  mainly  revolved  around  eompleting  the  implementation  of  support 
for  2.5-D  data  sets.  Issues  fixed  were  the  proper  display  of  eoordinate  slices  and 
generation  of  field  lines. 

•  Correeted  bugs  that  eaused  erashes  in  the  animation  module  for  Heliospaee 

•  Fixes  to  get  proper  graphieal  display  when  using  multiple  overlay  planes 

•  Expanded  the  2D  plotting  graphies  to  ereate  vertieal  labeling  with  proper  plaeement 

•  Created  2D  plot  graphing  object 

•  Created  satellite  grid  objeet  and  user  interfaee 

•  updated  various  models  for  V3.0  support 

•  updated  the  radiation  models  whieh  involved  updating  the  radiation  belt  eode  and 
ereating  APEX,  CURES  and  NASA  modules 

•  updated  magfield  and  SAA  models 

•  developed  a  new  grid  module  and  GUI  that  supports  various  eoordinate  spaees 

•  updated  the  satellite  grid  and  GUI 

Efforts  were  finalized  on  the  support  updates  to  various  models  for  V3.0.  These  updates  along 
with  all  pertinent  information  and  eode  were  sent  to  AFRL. 

Heliospaee  Analysis  Tools:  The  Heliospaee  software  paekage  is  a  unique  scientifie  tool  that 
provides  enhaneed  graphieal  and  data  visualization  of  3 -Dimensional  (3D),  time -varying 
systems.  The  software  was  developed  by  Boston  College  personnel  as  a  fundamental  graphies 
module  for  the  AFGEOSpaee  spaee  environment  program. 

The  Heliospaee  paekage  is  unique  in  that  it  generates  3D  adaptive  grids  for  both  Cartesian  and 
spherieal  data  sets.  In  addition,  field  lines  ean  be  generated  from  veetor  data  and  then  rendered 
into  3D  grid  spaee.  Graphieal  options  inelude  eoordinate  sliees,  grid  points,  isosurfaees  and 
domains.  Onee  rendered  the  graphics  may  then  be  manipulated  by  rotating,  translating  and 
zooming  in  and  out.  Multiple  viewports  ean  be  ereated  to  give  different  views  of  a  single  objeet 
or  to  view  multiple  data  sets  simultaneously.  In  addition,  time  dependent  data  sets  ean  be  viewed 
in  animation. 

During  the  eourse  of  this  eontraet,  the  Heliospaee  package  was  enhaneed  to  allow  a  user  to  ereate 
and  playbaek  maeros;  to  ereate  a  eolor  map  legend;  to  ereate  a  new  graphieal  objeet  to  allow  the 
user  to  add  annotations  to  the  image;  to  improve  display  options  for  the  data  domain;  and  to 
allow  the  user  to  seleet  the  baekground  eolor  and  gradient. 

3,4  History  of  the  Ring  Current 

The  "ring  eurrent"  grows  in  the  inner  magnetosphere  during  magnetie  storms  and  eontributes 
signifieantly  to  oharaeteristie  perturbations  to  the  Earth's  field  observed  at  low-latitudes.  Work 
under  this  eontraet  was  expended  on  writing  a  short  history  of  the  ring  eurrent  that  develops  in 
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the  Earth’s  magnetosphere  during  magnetic  storms.  This  effort  resulted  in  a  publication  in  the 
History  of  Geomagnetic  and  Space  Science  outlining  how  understanding  of  the  ring  current 
evolved  during  the  half-century  intervals  before  and  after  humans  gained  direct  access  to  space. 
Its  existence  was  first  postulated  in  1910  by  Carl  Stormer  to  explain  the  locations  and 
equatorward  migrations  of  aurorae  under  stormtime  conditions.  In  1917  Adolf  Schmidt  applied 
Stormer's  ring-current  hypothesis  to  explain  the  observed  negative  perturbations  in  the  Earth's 
magnetic  field.  More  than  another  decade  would  pass  before  Sydney  Chapman  and  Vicenzo 
Eerraro  argued  for  its  necessity  to  explain  magnetic  signatures  observed  during  the  main  phases 
of  storms.  Both  the  Stormer  and  Chapman-Eerraro  models  had  difficulties  explaining  how  solar 
particles  entered  and  propagated  in  the  magnetosphere  to  form  the  ring  current.  During  the  early 
1950s  Hannes  Alfven  correctly  argued  that  the  ring  current  was  a  collective  plasma  effect,  but 
failed  to  explain  particle  entry.  The  discovery  of  a  weak  but  persistent  interplanetary  magnetic 
field  embedded  in  a  continuous  solar  wind  provided  James  Dungey  with  sufficient  evidence  to 
devise  the  magnetic  merging-reconnection  model  now  regarded  as  the  basis  for  understanding 
magnetospheric  and  auroral  activity.  In  the  mid-1960s  Eouis  Erank  showed  that  ions  in  the  newly 
discovered  plasma  sheet  had  the  energy  spectral  characteristics  needed  to  explain  the  ring 
current's  origin  [36]. 

3,5  Second  Planar  Langmuir  Probe 

The  objective  of  this  effort  was  to  provide  engineering  services  for  a  number  of  circuit  boards 
and  other  elements  for  space  measurements.  One  of  our  efforts  included  investigating  the 
development  of  the  second  planar  Eangmuir  probe  (SPEP).  This  included  refining  the  design  of 
the  logarithmic  amplifier  and  laying  out  the  printed  circuit  board  for  it.  The  board  was  mostly 
laid  out  when  a  shortcoming  of  the  operational  amplifier  to  be  used  was  discovered.  The 
amplifier  input  has  to  be  limited  to  less  than  2.5  volts.  The  problem  was  discovered  in  testing  at 
Kirtland.  This  necessitated  much  circuit  redesign.  The  design  was  coordinated  with  Kirtland 
personnel. 

Pinal  schematics  of  the  analog  board,  board  layout  and  parts  list  for  the  engineering  version  of 
the  board  needed  were  prepared  for  presentation  to  APRL.  The  engineering  printed  circuit  board 
was  ordered  and  assembled  by  APRL  and  sent  to  Boston  College  for  testing.  Between 
engineering  board  design  and  testing,  time  was  spent  on  the  flight  version  of  the  board.  The 
flight  version  differs  in  the  component  footprints.  Some  commercial  components  were  used  for 
the  engineering  version  of  the  analog  board.  They  cost  a  small  fraction  of  the  flight  quality 
components.  Punctionally  they  are  the  same  but  have  different  packages  requiring  different 
printed  circuit  board  footprints.  Documentation  and  shipping  of  test  boards  was  completed. 

Work  was  also  done  in  collaboration  with  APRL,  MIT  Lincoln  Lab  (MITEL)  and  COSMIAC 
Corporation  in  the  design  of  a  new  radiation  detector.  The  new  design  instrument  was  to  be  a 
small,  light  weight,  low  power  radiation  detector.  Integration  details  between  APRL  hardware 
and  the  MITEL  satellite  were  completed.  The  interface  rules  kept  changing  as  MITEL  changed 
their  satellite  design. 

The  new  project  CEASE  Mk3  (Compact  Environmental  Anomaly  Sensor  Mark  3)  was  also 
addressed  under  this  contract.  The  Mk3  indicates  the  third  generation  of  similar  instruments.  The 
goal  of  this  instrument  is  to  provide  a  simple,  economical,  future  parts  available,  in-house 
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version  of  the  previous  generations.  The  previous  versions  had  many  unique  and  expensive 
parts.  Advances  in  the  electronics  will  permit  use  of  more  common  off  the  shelf  parts.  A  new 
power  supply  was  designed. 

3,6  Data  Analysis  Center  (DAC)  Support 

While  at  Hanscom  AFB,  Boston  College  team  members  supported  the  “Middle  Earth”  computer 
cluster  and  the  Stoermer  data  server  for  RVBX.  All  servers  utilized  Red  Hat  and  CentOS 
operating  systems.  Workstations  were  a  mix  of  Fedora  Linux  and  Microsoft  Windows  operating 
systems.  Team  members  performed  necessary  actions  in  compliance  with  all  security  and 
network  regulations.  System  performance  management  was  an  active  and  ongoing  task.  The 
team  also  maintained  hardware,  software,  and  software  licensing.  Research  of  equipment, 
peripherals,  and  software  was  included  in  the  overall  consulting  effort  and  recommendations 
were  provided  to  management  for  consideration.  Lastly,  user  support  demanded  a  significant 
portion  of  the  overall  effort. 

Throughout  the  BRAC  transitional  period,  Boston  College  provided  primary  support  at  AFRL 
Hanscom  AFB  for  archiving,  cataloging,  labeling,  and  properly  packing  servers  and  data  storage 
equipment  for  the  physical  move  to  Kirtland  AFB.  Upon  arrival  at  Kirtland,  BC  team  members 
provided  ongoing  technical  (interchange)  guidance  in  support  of  the  successful  reassembly  of  all 
RVBX  servers  and  data  storage  equipment  into  the  new  data  analysis  facility.  All  network 
communications  between  servers  and  storage  equipment  was  reestablished.  Operating  systems, 
security  and  network  settings  were  adjusted  for  the  new  base  guidelines  at  Kirtland,  and  all  data 
was  again  archived  and  verified.  BC  team  members  worked  tirelessly  through  all  challenges 
presented  by  the  BRAC  transition  and  until  all  tasks  were  successfully  completed. 

4,  RESULTS  AND  DISCUSSION 

Since  this  project  consisted  of  numerous  research  efforts,  we  have  included  the  results  and 
discussion  for  each  project  in  the  previous  sections. 

5.  CONCLUSIONS 

Under  this  contract,  the  Boston  College  ISR  conducted  an  innovative  scientific  research  and 
analysis  program  to  understand  the  fundamental  solar  and  magnetospheric  processes  and 
mechanisms  that  drive  and  control  the  Earth’s  space  environment  in  order  to  understand,  specify, 
anticipate  and  mitigate  the  effects  of  space  weather  on  Air  Force  Systems.  As  part  of  our  efforts, 
we  supported  this  program  with  research  and  instrument  development  to  study  the  origins  and 
interplanetary  (IP)  propagation  of  disturbances  that  affect  the  geosphere  and  supported  the 
mission  and  data  analysis  for  the  AFRL  Solar  Mass  Ejection  Imager  (SMEI)  experiment,  a 
Heliospheric  image/detector  that  observes  and  tracks  coronal  mass  ejections  from  the  Sun  to 
Earth.  Staff  members  also  contributed  to  the  development  of  the  Space  Weather  Forecasting  Lab 
(SWFL)  at  AFRL.  Boston  College  also  provided  scientific  and  technical  collaboration  to 
AFRL/RVB  to  aid  and  enable  the  development  and  validation  of  the  AE-9  and  AP-9  models,  the 
next-generation  of  trapped  electron  and  proton  models  designed  to  replace  the  current  AE-8  and 
AP-8  models  and  provide  users  with  much-needed  improvements  in  accuracy  and  additional 
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capabilities.  This  work  also  addressed  the  development  of  enhaneed  versions  of  AFGEOSpaee 
and  the  Heliospaee  analysis  toolset.  Finally,  special  projects  in  support  of  Earth’s  Space 
environmental  studies  were  addressed  throughout  the  duration  of  the  contract.  These  inelude  the 
use  and  analysis  of  various  space  sensors  and  satellite  missions,  and  the  development  of 
Stoermer  as  the  primary  data  file  server  for  RVBX.  Numerous  presentations  were  made  at 
national  and  international  conferenees.  At  least  20  papers  were  published  in  peer  reviewed 
journals. 
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AEB  Air  Eoree  Base 

AEGEOSpaee  Air  Eoree  GEOSpaee  Model 
AERE  Air  Eoree  Researeh  Eaboratory 

AESPC  Air  Eoree  Spaee  Command 

AEWA  Air  Eoree  Weather  Ageney 

AIA  Atmospherie  Imaging  Assembly 

AE-8  Trapped  Eleetron  Model  -  version  8 

AE-9  Trapped  Eleetron  Model  -  version  9 

AP-8  Trapped  Proton  Model  -  version  8 

AP-9  Trapped  Proton  Model  -  version  9 

APEX  Aetive  Plasma  Experiment  Satellite 

AU  Astronomieal  Unit 

BC  Boston  College 

BRAC  Base  Realignment  and  Closure 

CEASE  Compaet  Environmental  Anomaly  Sensor 

CD  AW  Coordinated  Data  Analysis  Workshop 

CME  Coronal  Mass  Ejeetion 

COSMIAC  Configurable  Spaee  Mierosystems  Innovations  and  Applieations  Center 

CURES  Coordinated  Radiation  Release  Satellite 

DAC  Data  Analysis  Center 
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DSX  Deployable  Struetures  Experiment 

EPAM  Eleetron,  Proton  and  Alpha  Monitor 

EUV  Extreme  Ultraviolet 

EVE  Extreme  Ultraviolet  Variability  Experiment 

EITS  Elexible  Image  Transport 

GOES  Geosynehronous  Operational  Environmental  Satellite 

GUI  Graphieal  User  Interfaee 

HAE  Hakamada-Akasofu-Ery  solar  wind  model 

HEO  High  Earth  Orbiter 

HE  High  Erequeney 

HI(s)  Heliospherie  Imager(s) 
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ID  IDentification 

IP  Interplanetary  Propagation 

IPS  Imaging  Proton  Sensor 

ISR  Institute  for  Seientifie  Researeh 

ISSI  International  Soeiety  for  Seientometries  and  Infometrics 
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LANL 
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Lm 
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MITLL 

Mk 
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NM 
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NRT 
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PAD 

PCA 

PI 

PRF 

RBSP 
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SDO 
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SEP 

SEPEM 

SEPMOD 

SKp 

SMEI 

SNRTACS 

SOHO 
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SPLP 

STEREO 

STOA 

STP 

SWEE 

TacSat4 

TB 

T-H 

TOA 

TSX5 

TWINS 

UCSD 
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Mega  Electron  Volt 

Multiple  Instrument  Chamber 
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Space  Dynamics  Observatory 
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Solar  Energetic  Particle 

Solar  Energetic  Particle  Environment  Modelling 
Solar  Energetic  Particle  Model 
Daily  average  solar  wind  speed 
Solar  Mass  Ejection  Imager 

Space  Nuclear  Radiation  Threat  Assessment  Code  System 

Solar  and  Heliospheric  Observatory 

Solar  Particle  Events 

Second  Planar  Eangmuir  Probe 

Solar  Terrestrial  Relations  Observatory 

Shock  Time  of  Arrival  Model 

Space  Test  Program 

Space  Weather  Eorecast  Eaboratory 

Tactical  Satellite  4 

Terabyte 

Tappen-Howard  Model 
Time  of  Arrival 

Tri-Service-Experiments  Mission  5  Satellite 
Two  Wide-angle  Imaging  Neutral-atom  Spectrometers 
University  of  California,  San  Diego 
Ultraviolet  radiation 


38 

Approved  for  public  release;  distribution  is  unlimited. 


Vh 

WHI 

WSA 

XML 


horizontal  plasma  drift 
Whole  Heliospheric  Interval 
Wang-Sheeley-Arge 
Extensible  Markup  Language 
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